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First principles density functional calculations for Fe2 03 has been performed over a wide range of 
pressures. The ground state is corundum-type hematite and is an antiferromagnetic insulator. This 
is in good agreement with experiment and other theoretical studies. On increasing pressure, the 
ground-state high-spin magnetic phase transforms to low spin via the closure of the charge transfer 
gap. The system also evolves to a new orthorhombic structure. Distorted corundum or Rh203(II) 
type structure with Pbcn symmetry and Pbnm type perovskite structure are two known competitive 
candidates for this structural phase, based on sm(/Ze-cationic type and tw o-cationic type picture, 
respectively. In our calculations, at about 38 GPa, Rh203(Il) type structure becomes more stable 
with respect to the ground state hematite. Relative stability of Pbnm type perovskite is ruled out 
by our calculations in this pressure regime. The Rh2 03(II) type structure remains in its low spin 
state, with 1 /xs/Fe atom, up to about 120 GPa. At this pressure the nonmagnetic solution in 
Rh203(ll) type structure becomes more favorable with respect to the low spin one. By further 
increasing the pressure at about 330 GPa, the system evolves to yet another new structural phase. 
This new orthorhombic structural phase is nonmagnetic and has Pmc2i symmetry, a subgroup of 
Cmcm. Surprisingly, on furthering rising the pressure, a Pbnm type nonmagnetic solution becomes 
competitive with the Pmc2i type structure and finally becomes stable at about 880 GPa. 
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I. INTRODUCTION 

Iron (Fc), one of the most abundant element on earth, 
is believed to contribute significantly to the mantle core 
characteristics in pure oxide form (FeO, Fe203) or in 
mixed oxide forms such as (Mg,Fe)0 and (Mg,Fe)Si03. 
Of these oxides, Fe203 shows versatility in it's structural, 
magnetic and electronic properties. The ambient pres- 
sure phase a-Fe203 (hematite) is corundum- type and 
is a wide-gap antiferromagnetic insulator^. It becomes 
weakly ferromagnetic between Morin temperature, Tm 
(=260 K), and Neel temperature, T^v (= 955 K), as a 
result of the canting of the spins of the two sublattices. 

Owing to its diverse intriguing properties it drew a 
significant amount of attention in the scientific commu- 
nity from the sixties. Reid and Ringwood^ based on Mc- 
Queen and Marsh's shock-wave experiment'^ proposed a 
new (denser) structural phase to exist between 60 - 120 
GPa. In the subsequent years the onset of an high pres- 
sure phase at 50 GPa was confirmed by quite a few 
studies'*''^'^'^. However, there was ambiguity concerning 
the crystal structure of the new phase. While Shan- 
non et. opted for the distorted corundum (Rh2 03- 
II) structure speculating on a single cationic nature of 
Fe, i.e. Ye^+ , the group of Suzuki et. al.^ went with the 
idea of orthorhombic perovskite structure with ABO3 for- 
mula with a ttiJo-type of cationic Fe (Fe^+ and Fe^+ ) pic- 
ture. Room temperature Mossbauer spectroscopy (MS) 
experiments'^'®'^ revealed a non magnetic component at 
~ 50 GPa, coexisting with the magnetic phase. A volume 
collapse of ~ 10%, concurrent with the onset of the phase 
transition at ~ 50 GPa, corroborating a first order phase 
transition, was observed by Olsen et. al. ' . Assuming the 
two cationic type picture this volume collapse was quali- 



tatively explained by the increase in coordination number 
(from 6 to 8) in one of the Fc sites. 

To resolve the issue of one or too-cationic types and 
correlate/separate the spin and structural transitions a 
number of investigations have been done in the more re- 
cent years. Of them, X-ray diffraction (XRD), MS at 
300K and electrical resistivity measurements by Paster- 
nak et .al.^^ unequivocally assigned high pressure struc- 
tural phase as the nonmagnetic distorted corundum or 
Rh203(II) type. They explained the intermediate region 
between the insulating and metallic region in their elec- 
trical resistivity data qualitatively by the coexistence of 
the insulating corundum phase (denoted HPl from now 
on) with the metallic Rh2 03(11) (denoted IIP2 in the 
following) . 

Other XRD and X-ray emission spectroscopy (XES) 
experiments by Badro et .al.^^ could separate the elec- 
tronic transition from the crystallographic transition. 
From their XES data it was concluded that the high spin 
(HS) could be stabifized in the high pressure structural 
phase, i.e. in IIP2 phase, at low temperatures. At high 
temperatures, however, the low spin (LS) state of IIP2 is 
stable. That this LS state is a weakly paramagnetic one 
was concluded by the presence of a satellite (K/3') in the 
XES spectra." 

The room temperature Raman spectroscopy study by 
Shim and Duffy^^ found the upper boundary for this 
phase transition (where no HPl signal is left) to be be- 
tween 54 - 56 GPa on compression. While decompress- 
ing, however, some (frequency) modes of the HP2 phase 
survive down to 25 GPa. The combined experimental and 
theoretical study by Rozenberg et .a/.^'^ also reported of 
the sluggishness in the transition from HPl to HP2. The 
onset of high pressure phase start at about ~ 40 GPa 
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and at 76 GPa they could find only the components of 
the high pressure phase. Their theoretical calculations at 
76 GPa predict the stability of the Rh203(II) [i.e. HP2) 
phase with respect to the orthorhombic perovskite phase 
(opv in the following). 

The XRD study by Ono et .al.^^'^''^ suggested that high 
temperature heating has a major role to play in control- 
ling the phase transition pressure. The phase transition 
pressure of ^ 45 - 50 GPa at room temperature drops 
down by about 15-20 GPa when the samples are heated 
to overcome the kinetic barriers hindering the transition. 
Their estimated transition pressure was ~ 26.2 GPa at 
300K. However, in this study they were not conclusive 
about the crystal structure of the high pressure phase. 

There are a few theoretical studies of the Fe2 03 sys- 
tem in the hematite phase. ^^'^^'^'^'^^ To the best of our 
knowledge, however, there are no theoretical studies con- 
cerning the low and high pressures phases and their sta- 
bilities. Neither the issue of transition from HS to LS is 
addressed. In this work we will shed light on these issues 
of the Fe2 03 system from first principles approach, based 
on density functional theory (DPT). 

The paper is organized as follows: after a brief discus- 
sion of the calculational details we will proceed to the 
next sections dealing with the results of our calculations. 
We will start with the ambient pressure phase. a-Fe2 03 
at the GGA level. After that we will proceed with two 
(orthorhombic perovskite, Pbnm type and Rh203(II), 
Pbcn type) competing candidate structures for the high- 
pressure phase following the hematite one. The next 
subsection will focus on the theoretical search for the 
experimentally found, but structurally unresolved, post- 
perovskite type phase. At the end wc will comment on 
the nature of the spin transition and the evolution of the 
structure. 



II. CALCULATIONAL DETAILS 

We have performed first-principles density functional 
calculations to analyze the spin and structural tran- 
sitions in FC2O3. Calculations were performed using 
the PBE generalized gradient approximation (GGA)'^". 
Variable-cell-shape molecular dynamics^^ was used for 
the full structural optimization at arbitrary pressures. 
Eight electrons (3d^4s^) per Fe atom and six electrons 
(2s^2p'*) per O atom were taken in the valence in the 
generation of the ultrasoft pseudopotential within the 
scheme of Rabe Rappe Kaxiras Joannopoulos ultrasoft 
pseudopotential^'^ . Total energy convergence of the wave 
function cut-off and kinetic energy cut-off (for the charge 
density) were checked. And the final values used were 
35 and 380 Ryd, respectively. The Brillouin zone inte- 
gration was performed using a Monkhorst and Pack^'^ 
mesh of 6x6x6 points. For the HPl phase (corundum 
structure) we used the rhombohedral cell consisting of 2 
formula units (f.u.) of Fe2 03. For the HP2 phase, and 
possible HP3/HP4 phase, orthorhombic cell consisting of 



4 f.u. of FC2O3 were used. 



III. RESULTS AND DISCUSSIONS 
A. hematite 

In the elemental rhombohedral cell, reported in Fig. 
1, four Fe atoms, spread along the z-axis, are grouped 
into two sets: two Fe (B-C) atoms at a shorter distance, 
denoted by {Fe-Fe)short and a pair (A-B or C-D) at a 
larger distance, denoted by (Fe-Fe)/ong- Depending on 
the orientations of the magnetization of the different Fe 
atoms, there can be four possible magnetic configurations 

(AFl = ABCD :: + + , AF2 = ABCD :: + - 

AF3 = ABCD :: + - -+ and FM = ABCD :: + + ++, 
where -I- and — correspond to opposite magnetization di- 
rections). The optimized lattice parameters as extracted 
from energy-volume curves are 10.277, 10.317, 10.340, 
10.355 a.u. respectively for the AFl, AF2, AF3 and FM 
configurations. The most stable magnetic configuration 
is AFl, as can be seen from the curves in Fig. 1. These 
lattice parameters were then used as the starting point 
for the subsequent cell optimization at finite pressures. 
In all antiferromagnetic configurations the Fe magnetic 
moments are close to 3.8 jU^/atom at GPa. Moreover 
all antiferromagnetic solutions are insulating while the 
ferromagnetic solution is metallic. 

The partial density of states at GPa, for AFl, is 
shown in the lowest panel of Fig. 2. From the figure 
it is evident that the top of the valence band is occu- 
pied by the 2p states of O and the 3d states of Fe on 
an equal basis. The valence and conduction band edges 
are seperated by a well defined gap of about 0.3 eV. It 
is to be noted that, all the antiferromagnetic configura- 
tions show a similar valence and conduction edges, with 
similar value of the fundamental gap, while the FM con- 
figuration is metallic. The experimental scenario is a bit 
different: the top of the valence band has mainly 0-2p 
character, and the experimental value of the band gap 
is 2.0 eV [24]. Both these deficiencies can be attributed 
to an incomplete treatment of electronic correlation at 
the LDA or GGA level. Similar results have been re- 
ported earlier and, for instance, the calculated values of 
the band gap is 0.3 eV by employing pseudopotential^*', 
0.75 cV by augmented spherical wave method^'^ and 0.51 
eV by linear muffin-tin orbitaP^, respectively. The band 
gap value in a Hartree Fock calculation^^ is instead of 
grossly overestimated as 11 eV due to the complete ne- 
glect of correlation effects. A good representation of the 
band gap and of the atomic decomposition of the band 
edge states can be obtained by including strong correla- 
tion effects via GGA + U, as can be seen in the upper 
panels of Fig. 2. We wish to emphasize here, however, 
the qualitative agreement of the GGA results with ex- 
periments of the basic electronic structure of hematite. 

To check the relative stability of the different magnetic 
configurations af finite pressures wc have calculated their 
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FIG. 1: (color online) Left panel: The rhombohedral unit 
cell of Fe2 03. The (big) green balls correspond to Fe 
atoms and brown ones correspond to Oxygen atoms. Right 
panel: Energy- Volume curves (symbols correspond to calcu- 
lated data and the curves are fits to data obtained using Mur- 
naghan's equation of state) for the four magnetic configura- 
tions of hematite considered in this work. Note that the curve 
corresponding to AFl configuration (see text for notation) has 
the lowest energy at the minimum and is therefore the ground 
state configuration. 
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FIG. 2: (color online) Projected Density of States for the 
hematite (AFf ) at GPa per atom for GGA and GGA -|~ U 
calculations. 



enthalpies by fully optimizing each of the structures by 
variable-cell damped dynamics at several pressures. The 
results are collected in Figure. 3 



B. spin transition 

As can be seen in Fig. 3(a), for all the magnetic phases 
the value of the magnetic moment per Fe atom remains 
almost constant at about 3.8 ^s/Fe atom at low pres- 
sure (large volumes). With decreasing volume (increasing 
pressure) there is a sudden collapse in magnetic moment 
magnitude to ^ Ifis/Fc atom. 

The relative enthalpy curves of HS and LS phases for 
the various configurations are depicted in Fig. 3(b). The 
crossover from the HS to LS occurs at around 36, 24.5 
and 19.5 GPa, for the AFl, AF2, AF3 configurations, re- 
spectively, and at 5 GPa for the FM (not shown in the 
figure). From this figure it is clear that the AFl config- 
uration is the most stable one upto a pressure of about 
34 GPa, where the FM phase becomes more stable. This 
trend was also observed by RoUman et.al..^'' Overall, the 
computed sequence of relative stability of the different 
magnetic phases of hematite is in agreement with exper- 



iment and previous calculations. ^^'^^ 

This HS-LS spin collapse at the transition is also ac- 
companied by a sudden change in the magnitude of a 
number of inter- and intra- octahedral structural param- 
eters. Notable inter-octahedral structural parameters are 
the {Fe~Fe)short and (Fe-Fe)/orag distances, which de- 
crease and increase, respectively, at the spin transition 
as can be seen from Table I. The same was observed by 
RoUman et.al.^^ in their GGA calculations. 

Although also the remote surroundings {i.e. the next 
to next-nearest neighbors and inter-octahedral distances) 
have non-negligible effects on the nature of the spin tran- 
sition, the most significant changes observed at the tran- 
sition concern mainly the octahedra, as already pointed 
out by a number of studies. Of the intra- 
octahedral structural parameters, the average Fe-0 bond 
lengths, as shown in Fig. 3(c), show a jump from higher 
to lower values and display two very distinct and well 
separated slopes in the HS and in the LS configurations. 

Two other intra-octahedral structural parameters that 
displays significant changes at the transition are the Fe- 
O bond-length distortion and the octahedral-angle vari- 
ance. Fe-0 bond-length distortion, Sdpe-Oi is defined 
1 ^ 

as Sdpe-o ~ — l^o'ilj where D is the average (octa- 

i=l 
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TABLE I: Variation of Fe-Fe bond lengths and of the average Fe-O bond lengths in the octahedral coordination across the spin 
transition. Increases (j) and decreases(|) in bond lengths are marked by arrows. The values of the distortion index, Sdpe-o, 
and of the octahedral angle variance, A9, at the transition are also reported. 
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FIG. 3: (color online) Results for the different magnetic con- 
figurations in the hematite structure, (a) Magnetic moment 
as a function of volume calculated at the GGA level, (b) 
Relative enthalpy curves. Solid lines (vertical) correspond to 
high spin (HS) and dashed lines (horizontal) are for low (LS: 
1/is/Fe atom) spins. The reference enthalpy (base line) cor- 
responds to the FM configuration in the LS phase, (c) Depen- 
dence of the average Fe-O bond-lenght (in the octahedra) with 
pressure. Vertical dotted lines correspond to the spin transi- 
tion pressures for different magnetic configurations. For the 
FM configuration, not shown in figure, similar characteristics 
are observed but at much smaller pressure of around 5 GPa. 



hedral) Fe-O bond-length at that particular volume and 
Adi is the difference of the i-th Fe-O bond-length from 
its average. The octahedral-angle variance''^ is instead 

12 

given by M'^ ^ ^{^j - 90)^/11. Where Bj is the vari- 

ation of the j-th 0-Fe~0 angle from its ideal value of 
90°. 

At the spin transition, both the mentioned intra- 
octahcdral structural parameters, Sdpe-o a-nd A0^, show 
a jump from a larger to a smaller value as can be seen 
from Table. I. This indicates that at the spin transition 
the octahedra become more regular and this happens ir- 
respective of the octahedral surroundings.'^'^ Moreover, 
the sudden jump of these parameters can directly be re- 
lated with the ~10% volume collapse that is observed for 
all the magnetic phases at the spin transition. 



C. structural transition 

Along with the spin transition discussed above, a- 
FC2O3 also undergoes structural phase transition/s un- 
der pressure and temperature. Two candidate structures 
for this phase are the orthorhombic perovskite {opv) with 
Pbnm symmetry and the Rh203(II) with Pbcn symmetry. 

In hematite (a-Fe203), in the corundum structure, the 
Fe atoms sit inside the cage of a single kind of poly- 
hcdra which are octahedra. Each of these octahedra 
shares three edges and one face with the neighbouring 
octahedra. Rh2 03(11) structure with Pbcn symmetry is 
of the distorted corundum type and also has the same 
single type octahedral arrangement. Owing to a larger 
octahedral distortion and rotation of the octahedra in 
Rh203(II), one of the three sharing edges is broken by 
the alteration of one of the coordinating oxygens. This 
crystallographic change is minor as can be seen from the 
Fig. 4: instead of three shared edges as in hematite, oc- 
tahedra in Rh203(II) have two sharing edges with the 
addition of a corner sharing bonding. 






FIG. 4; (color online) Polyhedral structures of Fe203 in dif- 
ferent crystal packing, (a): hematite, (b): Rh203(II), (c): 
Pbnm type perovskite. Small white spheres are oxygen atoms 
and big green ones are for the Fe atoms. 

On the other hand, for the pervoskite structure with 
Pbnm symmetry, there are two types of polyhedral en- 
vironment surrounding Fe atoms. The six-fold coordi- 
nated (octahedra) and the 8-fold coordinated polyhcdra 
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are arranged in alternate layers along the z direction. 
The octahedral connectivity is via corner sharing, while 
the 8-fold coordinated polyhedra arc connected to each 
other by edge and corner sharing. The interpolyhedral 
connectivity (between octahedra and 8-fold coordinated 
polyhedra) is maintained via face and edge sharing. 

A Mossbaucr spectroscopy study at 300K by Paster- 
nak et. al}'^ showed that the transition from hematite 
to the high-pressure phase starts at about 45 GPa and 
completes at about 72 GPa. The existence of o-pv phase 
at high-pressure was ruled out by pointing out the single 
cationic nature of the Fe-ions. Moreover, Fe-ions were 
characterized as non-magnetic and the magnetic to non- 
magnetic transition connected to the structural phase 
transition. In a later study Badro et. al}^ could separate 
the crystallographic transition from the spin transition. 
According to this study the spin transition is isostruc- 
tural and takes place in the high-pressure phase. In a 
more recent study Ono et. al}^'^'^ reported the transition 
pressure to be 26.2 GPa at 300 K. It was concluded there 
that high-temperature heating played a fundamental role 
in overcoming the kinetic barrier hindering the phase 
transition, thus lowering the transition pressure. The 
combined experimental and theoretical study by Rozen- 
berg et. al}'' also indicated Rh2 03(11) as the IIP2 phase. 
The o'pv phase was found to be higher in energy than the 
Rh203(II) phase by as much as 25 Ryd/f.u. at 76 GPa. 
However, no other detailed theoretical study exist on the 
relative stability of these two phases. 

We have calculated the static enthalpies for these two 
structures, both in FM and AF configurations, for a finite 
range of pressures. Fig. 5 shows the relative enthalpy- 
pressure curves for these two candidate structures along 
with the AFl corundum structure. As can be seen from 
the figure, at low pressures all structures show an HS 
solution. Corundum in the HS AFl configuration is the 
most stable one upto about 34 GPa. At this pressure 
the enthalpy of corundum in the LS FM configuration 
becomes lower. This crossover from HS AF to LS FM is 
also observed for Rh2 03(11) structure at about 26 GPa 
(as can be seen from the crossing of the HP2-af-hs and 
HP2-fm-ls curves in Fig. 5). The crossover from HS to 
LS is also present for opv structure, but only in the oc- 
tahedral sites. The 8-fold coordinated polyhedra remain 
in the high spin and hence, overall, the system remain 
in a kind of intermediate mixed (half HS and half LS) 
spin state. We mention also that from ~ 34 GPa to 
about 38 GPa the LS FM configuration in the corun- 
dum structure is the most stable one. After this, in the 
broad range of pressure from ~ 38 GPa to almost 120 
GPa the Rh2 03(11) structure in the LS FM configura- 
tion is the moste stable. With increasing presssure there 
is a gradual decrease of magnetization in the system that 
finally becomes non-magnetic from ^ 120 GPa onwards. 
Hence, upto the pressure range shown in this figure, opv 
does not have any stability region. Therefore, our cal- 
culations, in agreement with recent experiments, ^'^'^^'^'^ 
rule out the two cationic site picture, i.e. opv, and estab- 



lishes Rh203(II) as the HP2 or next stable phase after 
the corundum one.. 
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FIG. 5: (color online) Relative enthalpy curves for the 
Rh2 03(11) (HP2), opv (opv) along with the corundum struc- 
ture (HPl). HPl structure in the LS phase is dynamically 
stable upto about 180 GPa. For the opv structure in the FM 
configuration, not shown here, the relative enthalpy curve re- 
main well above the base line here. The dashed vertical lines 
are guide to the eye to locate the transition pressures. 



D. further structural transitions 

The high temperature experiments by Ono et. al.^'^ 
reported of another structural phase transition after the 
HP2 phase. Their data show that the transformation of 
this phase begins at around 60 GPa and 2500K. They 
identified this phase as Calr03 type with Cmcm sym- 
metry. However, this designation was made without a 
detailed structural analysis in analogy to other systems. 
The prototype scsquioxidc, AI2O3, for instance, having 
corundum structure at ambient conditions, undergoes a 
structural transformation from Rh2 03(11) type phase to 
a post-perovskite one of Calr03 type, and this structure 
has been found as the post-perovskite phase in most of 
the instances of post-perovskite phases discovered till to- 
day [see Ref.34 and references therein]. 

In FC2O3, being another sesquioxide, it is therefore 
natural to explore the possible higher pressure phases 
starting with Calr03 structure with Cmcm symmetry 
(denoted as HP3 from here). In terms of polyhedral ar- 
rangements this structure has two types of polyhedral ar- 
rangements, i.e. octahedra and 8-fold coordinated poly- 
hedra. The difference with respect to the Pbnm type 
structure is that there is no longer any inter-layer octa- 
hedral connection. Instead alternate layers of octahedra 
and 8-fold coordinated polyhedra form the structure as 
shown in Fig. 6(a). 



FIG. 6: (color online) Two candidate high pressure structures 
of Fe203. (a) HPS: Cmcm type post-perovskite and (b) HP4: 
Pmc2i type post-perovskite. Small white spheres are oxygen 
atoms and big green ones are for the Fe atoms. 
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In Fig. 7, one can see that around 200 GPa both opv 
and the HPS phases are in the intermediate spin state 
(Fe in the octahedral cage is in low spin state, while it 
is in a high-spin state when 8-fold coordinated to Oxy- 
gen) and are dynamically stable. However, upto about 
337 GPa, the enthalpy of the nonmagnetic IIP2 phase 
is lower. From 337 GPa onwards, the enthalpy of the 
nonmagnetic phase of a new structure''''' crosses the IIP2 
one. This new structure, denoted HP4 in the following, 
has Pmc2i symmetry with an orthorhombic unit cell con- 
sisting of 20 atoms. Pmc2i is a subgroup of Cmcm. The 
4c, 4a and 8f Wyckoff positions in Cmcm are split into 
2a and 2b Wyckoff positions in Pmc2i space group. This 
splitting of the symmetry positions leads to different kind 
of polyhedral arrangement in the new phase. Instead of 
alternate layers of octahedra and of 8-fold coordinated 
polyhedra as in Cmcm, the arrangement in Pmc2i cor- 
responds to alternate layers of octahedra and of 7-fold 
coordinated polyhedra. Compariing HP4 with IIP3, it 
can be seen from Fig. 7 that the AF IIP3 curve, labelled 
HP3-af-is, crosses the non-magnetic HP4 curve {HP4- 
non-mg) at ~ 225 GPa. This newly achieved structural 
phase (HP4) remains stable for a wide pressure range. 
At much higher pressures, however, the enthalpy of the 
nonmagnetic opv phase becomes competitive and finally 
becomes the stable structure at about 880 GPa. 

Experimentally,^'* a post-perovskite type phase was 
obtained at pressure ( 60 GPa at 2500K) much be- 
low the one achieved here by our calculations, where 
the post-perovskite phase is stable with respect to the 
Rh203(II) phase only beyond 330 GPa. High tempera- 
ture heating likely plays a role in reducing the experi- 
mental transition pressure. For instance, for the struc- 
tural transition from corundum- type to Rh2 03(11), the 
transition pressure drops from 46 GPa to 26 GPa when 
the samples are laser heated from room temperature to 
about 2500K.^*'*'''^''' We notice that the Cmcm type post- 
perovskite structure proposed (without structural deter- 
mination) by Ono. et. a/. is never the stable one in 



FIG. 7: (color online) Relative enthalpy curves for the 
Rh203(II) (IIP2), opv (opv), Cmcm type post-perovskite 
(HPS) and Pmc2i type post-perovskite (HP4). The dashed 
vertical lines are guide to the eye to locate the transition pres- 
sures. 



our calculations, rather Pmc2i symmetry-a subgroup of 
Cmcm-is found. Further efforts are needed to establish 
experimentally the structural parameters of this post- 
perovskite phase. 

This extreme condition phase might have implications 
in exploring and understanding the rocky cores of solar 
giants and/or rocky mantles of recently discovered ter- 
restrial exoplanets."'^'"'^''''^. 



E. in the vicinity of spin transition and structural 
evolution 

Before coming to conclusions we want to summarize 
two issues: What happens at the spin transition and 
what is the overall structural evolution in Fe2 03. We 
have already mentioned that the phenomena of spin tran- 
sition occurs in the octahedral cages''^ for all the three 
phases that arc dynamically stable at ambient pressure: 
i.e. hematite, Rh203(II) type Fe203 and Pbnm type 
Fe203. A vivid look into the density of states profile, 
depicted in Fig. 8, shows that at the transition, the d~ 
band width increases substantially. Our explanation for 
this goes as: with increasing pressure the O atoms gets 
closer to the Fe atoms. This induces enhancement in the 
octahedral crystal field splitting of the triply degenerate 
t2g and doubly degenerate eg levels. At some point in 
the pressure scale the strength of the crystal field energy 
overcomes the exchange energy and the Fe atoms transit 
from a state with '^t2g ^Cg configuration to a state with 
^t2g °eg. Also, it can be seen from Fig. 5 that in terms of 
magnetic state the system exists in three distinct regions, 
i.e., with HS, LS and nonmagnetic. The transition from 
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FIG. 8: (color online) Projected density of states in the vicinity of the spin transition. Lower panels show the DOS just before 
the spin transition occurs and the upper panels are the same just after the spin transition, (a) for the hematite structure, (b) 
for the Rh203(II) type Fe203. (c) for the Pbnm type Fe203. Note that in Pbnm type Fe203 Fe is in two types of polyhedral 
coordination. The spin transition occurs only in the octahedral cage. Whereas the Fe in 8-fold coordinated polyhedra remain 
in high spin. 




(a) (b) (c) 

FIG. 9; (color online) Figure shows the evolution of the basic network of the polyhedra. E and F correspond respectively 
to edge sharing and face sharing bonding. Red (smaller) spheres are O atoms while the bigger spheres at the centres of the 
polyhedrons correspond to Fe atoms and are denoted by 1, 2, 3, 4. (a) correspond to Hematite, (b) Rh203(II) structure, and 
(c) Pmc2i structure. Note the changes in the edge, face and corner sharing, leading to different types of structures. 
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LS to the nonmagnetic state is rather gradual. This hap- 
pens only when the hybridization strength between Fe-3d 
and 0-2p levels overcomes the crystal field energy. 

To explain how the structure evolves, we present the 
sequence of basic polyhedral network from hematite to 
Pmc2i type structure through Rh2 03(11) structure in 
Fig. 9. The face sharing network in between ocathedron 
1 and 2 in hematite is retained in Rh203(II). However, 
the edge sharing network between octahedron 2 and 3 of 
hematite is broken and instead a corner sharing bonding 
is formed between octahedron 3 and 4 in Rh2 03(11). In 
the transition from Rh203(II) to the Pmc2i structure, 
the face sharing network between octahedron 1 and 2 of 
Rh203(II) is broken and the rearrangement of the inter- 
nal parameters leads to two kinds of polyhedral connec- 
tivity, one with 6-fold coordination (2 and 4) and the 
other with 7-fold coordination (1 and 3). 

IV. CONCLUSION 

In this work, systematic density functional calcula- 
tions of Fe2 03 have been performed over a wide range 
of pressure. Corundum type phase displays high-spin to 



low-spin transition for all magnetic configurations. This 
corundum type phase transforms to a distorted corun- 
dum or Rh2 03(11) type phase at about 38 GPa while 
Pbnm type perovskite structure don't have any stability 
field in this pressure regime. The phenomenon of spin 
transition is also present in Rh2 03(11) structure. On in- 
creasing pressure, our calculations identify a new struc- 
ture with Pmc2i symmetry, a subgroup of Cmcm, as the 
post-perovskite structure, while the experimentally pro- 
posed Cmcm phase is never the most stable one. Our 
calculations also show that the orthorhombic perovskite 
type structure with Pbnm symmetry might have a sta- 
bility field at extremely high pressures. 
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